arpa  order  no.  lse-ei 


MEMORANDUM 

RM-3505-ARPA 

FEBRUARY  1963 


THE  PROBABILITY  DISTRIBUTION  OF 
ANTI-MISSILE  MISSILE  MISS  DISTANCE 

DUE  TO  OBSERVATION  AND 
GUIDANCE  NOISE 

W.  B.  Kendall 


PREPARED  FOR: 

ADVANCED  RESEARCH  PROJECTS  AGENCY 


% 


SANTA  MONICA  .  CALIFORNIA 


ARPA  ORDER  NO.  180 -61 


MEMORANDUM 

RM-3505-ARFA 

FEBRUARY  1063 


THE  PROBABILITY  DISTRIBUTION  OF 
ANTI-MISSILE  MISSILE  MISS  DISTANCE 
DUE  TO  OBSERVATION  AND 
GUIDANCE  NOISE 

W.  B.  Kendall 


This  research  is  supported  and  monitored  by  the  Advanced  Research  Projects  Agency 
under  Contract  No.  SD-79.  Any  views  or  conclusions  contained  in  this  Memorandum 
should  not  be  interpreted  as  representing  the  official  opinion  or  policy  of  ARPA. 
Permission  to  quote  ®*r  reproduce  must  be  obtained  from  The  RAND  Corporation. 


nil  n  D, 


17«0  MAIN  ST  •  SANTA  MONICA  ■  CALIFORNIA  « 


iii 


PREFACE 


The  study  presented  in  this  Memorandum  steins  from  a  problem 
encountered  in  setting  accuracy  specifications  for  the  observing 
instruments  of  a  midcourse  intercept  system  intended  to  operate  against 
submarine -launched  ballistic  missiles.  The  result  given  here  is  much 
more  general  than  the  problem  stated,  and  can  be  applied  to  a  wide 
variety  of  intercept  situations.  The  study  is  part  of  RAND' s  continuing 
work  on  the  interception  of  submarine -launched  and  air-launched 


ballistic  missiles. 


V 


SUMMARY 


In  this  Memorandum  we  consider  the  problem  of  intercepting  pn 
attacking  ballistic  missile  (BM)  with  an  anti-missile  missile  (AM’’) 
when  Gaussian  noise  corrupts  both  the  observations  of  the  BM  and  the 
guidance  system  of  the  /'MM. 

It  is  essumed  that  the  position  of  the  BM  at  some  future  time  is 
estimated  from  a  set  of  (possibly  correlated)  measurements  (of  any 
type )  made  on  it,  and  that  the  error  in  this  position  estimate  is 
related  to  the  errors  in  the  measurements  by  the  usual  linear  relation 


x. 

l 


e  . 

<j 


where  x.  is  the  error  in  the  i-tn  special  coordinate  w.,  and  e.  is 
i  i  j 

the  error  in  the  j-th  measurement  m..  Furthermore,  it  is  assumed 

J 

that  the  noise  in  the  AMM  guidance  system  is  such  that  the  uncertainty 
in  the  position  of  the  AMM  at  intercept  has  a  general  three-dimensional 
Gaussian  distribution.  For  these  general  assumptions  the  probability 
that  the  miss  distance  between  the  AMM  and  the  BM  will  be  less  than 
eny  given  amount  is  determined  and  shown  in  several  curves.  In  the 
course  of  this  calculation  the  cumulative  distribution  function  for 
the  generalised  chi- squared  distribution  in  eny  number  of  dimensions 
is  obtained. 

For  the  more  special  case  of  un correlated  measurements  and  AMM 
position  uncertainty  with  spherical  symmetry  (this  might  be  the  case 
when,  for  example,  only  the  rms  miss  distance  of  the  AMM  rather  than 
its  complete  three-dimensional  distribution  is  known)  coefficients 


<? 


vi 


ere  obtained  which  indicate  how  sensitive  the  miss  distance  is  to  in¬ 
cremental  changes  in  the  accuracies  of  any  of  the  measurements.  With 
these  coefficients  a  system  designer  can  systematically  determine  how 
much  money,  weight,  and  equipment  to  allot  to  making  each  of  the 
measurements.  A  numerics!  example  is  given  which  demonstrates  this 


procedure. 
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1 

r 

c 
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i 
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s* 

cp l(v) 
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characteristic  function  of  f(r  ) 
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I.  INTRODUCTION 

Cons  id  erst  le  interest  has  been  shown  in  the  following  problem: 
A  set  of  measurements  is  made  on  an  attacking  ballistic  missile  (BK), 
and  based  on  these  measurements,  its  position  at  some  future  time  t 
is  estimated.  Due  to  Gaussian-distributed  errors  in  these  measure¬ 
ments,  the  error  in  the  position  estimate  has  a  Gaussian  distribution. 
An  anti-missile  missile  (AMM)  is  aimed  to  intercept  the  BM  at  time  t. 
However,  due  to  noise  in  its  guidance  system,  its  position  at  time 
t  is  in  error  by  a  random  amount  which  has  a  Gaussian  distribution. 
Determine  the  probability  that  the  distance  between  the  BM  and  the 
AMM  at  time  t  is  less  than  some  magnitude  k. 

In  this  Memorandum  we  attack  and  obtain  a  good  approximate  solu¬ 
tion  for  this  problem.  Of  even  more  importance,  in  the  solution  ob¬ 
tained  the  effects  of  AMM  guidance  noise,  of  the  tactical  geometry, 
and  of  the  measurement  accuracies,  can  be  separated.  Thus  we  nave 
a  simple  method  for  determining  the  change  in  kill  probability  due 
to  changes  in  any  of  these  variables. 
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II.  MATHEMATICAL  FORMULATION 


We  will  use  matrix  notation  throughout.  A  superscript  T  will 
denote  the  transpose. 

Let  the  position  estimate  of  the  BM  he  denoted  by  the  set  of 
rectangular  coordinates  -jw^j.  ~  v*  Each  of  these  is  in  general  8  func¬ 
tion  of  the  set  of  measurements  =  m.  If  the  have  a  jointly 

Gaussian  distribution  with  covariance  matrix  M,  then  to  what  is 
usually  a  good  approximation,  the  errors  in  w  (call  them  -jx^j.  =  x) 
are  jointly  Gaussian  with  covariance  matrix  X  =  DMD  .  Here  D  is 
the  "sensitivity  coefficient  matrix1'  given  by  D  =  {Dij}  v^ere 

dw . 

D .  .  -  -3—- 
1J 


Denote  by  =  z  the  guidance-noise  induced  errors  in  the  rec¬ 

tangular  coordinates  of  the  AMM  at  time  t.  These  errors  are  assumed  to 
to  be  jointly  Gaussian  random  variables  with  zero  means  and  covariance 
matrix  Z.  Now,  since  the  AMM  was  aimed  at  point  w,  its  position  at 
time  t  is  w  +  z.  The  actual  position  of  the  BM  at  time  t  is  w  -  x 
(i.e.,  the  estimated  position  equals  the  actual  position  plus  the  error). 


*The  i,j-th  element  of  the  covariance  matrix  is  the  expected  value 
of  eiej  where  ek  is  the  k-th  random  variable.  It  is  a  symmetric  matrix. 

This  is  because  it  is  usually  satisfactory  to  make  the  linear  ap- 

j.  is  the  matrix  of  measurement  errors  and 

dw. 

c  =  {Dij}  =  {k:}  is  the  sensitivity  coefficient  matrix.  Then  one  has 
J 

X  =  E(xxT)  =  E(DeeTDT)  =  D[E(eeT)]DT  =  DMDT  where  E  is  the  operator  which 
replaces  each  element  of  a  matrix  by  its  expected  value. 


proximation  x  -  De  where  e  -  J e 
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Thus  the  rectangular  coordinates  of  the  difference  in  their  positions 
is  x  +  z  =  y;  which  is  Gaussian  with  zero  mean  and,  assuming  statistical 
independence  of  w  and  z,  covariance  matrix  X  +  Z  =  Y. 

The  square  of  the  distance  between  the  BM  and  the  AMM  is 

r  =  £  (x±  +  zi)  =  y  y 
i 

T  2 

Our  problem  then  is  to  calculate  the  probability  that  y  y  <  k  for  any 
k,  given  that  the  probability  density  of  y  is 

•  f(y)  =  exp  j-  yT  Y"1  yj 

where 

=  (Zjr)"^2  (det  Y 

and  q  is  the  dimensionality  of  the  space  under  consideration.  Here  we 

O 

are  interested  in  q  =  2  or  3.  The  quantity  r  has  the  so-called  "gen¬ 
eralized  chi-squared"  distribution,  and  has  been  treated  by  many  authors, 

(2) 

most  notably  perhaps,  by  J.  I.  Marcum.  '  We  will  make  much  use  of 
thi s  work  of  Marcum T  s . 
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III. 


THE  CHARACTERISTIC  FUNCTION  OF  r 


2 


The  characteristic  function  of  r  , 


denoted  by  <P(v) , 


■* 

is  given  by 


O' 

cp(v)  =  J  exp  |jvr  (y)j  f(y)  dy 


=  J*  exp  |  jv  yTy  -  ^  /  Y"1  yj  dy 


=  (2n)-Ci//2  (det  Y)"1//2  J  exp  |  -  ~  yT  (Y”1  -  2jvl)  yj  dy 


det  (I  -  2 jvY ) 


q 

TT  (l  - 

i=l 


2jv  n_.) 


1/2 


where  the  are  the  reel,  positive  characteristic  roots  of  the  positive 
definite  covariance  matrix  Y.  They  are  the  solutions  of  the  equa¬ 
tion 

det  (Y  -  pi)  -  0 


where  "det”  denotes  the  determinant. 


j.  00  00 

By  p  dy  we  mean  j*  dy^  . . . 

-00  •“•00 


In  the  following  treat¬ 


ment  I  is  the  unit  matrix  which  has  ones  on  its  main  diagonal  and 
zeros  elsewhere. 
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IV,  THE  PROBABILITY  DENSITY  FUNCTION  OF  r2 

2 

The  probability  density  function  of  r  can  be  found  by  taking 
the  inverse  Fourier  transform  of  <p(v). 

f(r^ )  =  2~~  J*  e'JVr  cp(v)  dv 


=  —  r 

2jt  J 


.  2 
-jvr 


-«  C(1  ’  2jvn1)...(l  -  2jvp^}] 


1/2 


dv 


When  q  =  2  (the  plane  case)  this  can  be  integrated  in  closed  form  to  give 


f(r2)  =  \  (p1p2)"l/2  exp 


{-  i  (i  *  h) 


r  ^ 


-  +  -  i  r  ]  r  >  0 

hi  h2 


For  q  >  2  the  integral  cannot  be  evaluated  in  closed  form,  but  following 

2 

Marcum  we  can  expand  f(r  )  in  s  Laguerre  series,  the  coefficients  of 

2 

which  can  be  expressed  in  terms  of  the  moments  of  f(r  ).  These  can 
easily  be  determined  directly  from  <p(v)*  The  Laguerre  series  is 


f(r2)  =  £  c  e-t  ta  La(t) 
i=0 


t  >  0 


where 


t  =  r  /b 


La(t)  are  the  generalized  Laguerre  polynomials  and  a  and  b  are  arbitrary 
constants. 


For  the  integration,  see  Ref.  3>  pair  555* 
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By  multiplying  "both  sides  of  this  equation  by  L^(t),  integrating  from 

J 

zero  to  infinity,  and  using  the  orthogonality  property  of  the  Laguerre 

polynomials,  we  obtain  the  coefficients  c^  8S  functions  of  the  moments 
2 

of  f(r  )  and  the  arbitrary  constants  s  and  b.  By  choosing 


a  = 


v2  "  V1 


1 

-  1  =  -5-  -  1 

a 


and 


b  - 


V2  "  V1 


where  is  the  i-th  moment  about  the  origin  of  f(r  ),  we  get 


1  1 

C0  “  b  r(e  +  1)  2  /  2;2, 

CT  r(v,  /o  ) 


ci  =  c2  =  0 

and 


C3  b'r(a  +  4) 


(a  +  3) 


In  evaluating  these  coefficients  we  will  make  use  of  the  fact  thal 
the  moments  are  related  to  the  cuniulants  k^  by 


V1  “  kl 


v2  =  k2  +  kl 


V3  ~  k3  +  3klk2  +  kl~ 
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The  cumulants  axe  defined  "by 


Now,  since 


00  t 


<P  (jv)  k 
In  cp(v)  =  ^  — r~: - 

i=l 


cp(v)  =  det  (I  -  2jvY ) 


-1/2 


we  can  write 


In  cp(v)  =  -  --  In  det  (I  -  2jvY) 


.Next  we  make  use  of  the  identity 


det  (I  -  2jvY)  h  exp  j“  ^ 

i=l 


C24y.L  tr  Y1 


where 


to  obtain 


tri1  =£(,/ 


m=l 


*(v)  -\l 
i=l 


In 


(Here  tr  denotes  the  trace  operator.  Besides  this  representation, 
we  also  have  tr  Y  equals  the  sum  of  the  diagonal  elements  of  Y.) 


This  representation  is  valid  for  all  v  <  l/(2q  )  where  q  is 
the  largest  characteristic  root  of  Y.  ®  ® 
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Equating  coefficients  of  v1  in  the  two  expressions  for  In  cp(v)  we  see 
that  here 

\  =  21'1  (i  -  i):  tr  y1 

and  that 

Vf  =  tr  Y 

v2  -  2  tr  Y2  +  (tr  Y)2 

=  8  tr  Y3  +  6  tr  Y  tr  Y2  +  (tr  Y)3 

These  can  he  used  to  evaluate  the  coefficients  in  the  Laguerre  series. 

2 

However,  as  will  he  seen  below,  a  satisfactory  approximation  to  f(r  ) 
is  obtained  even  when  c^  and  the  higher  coefficients  are  neglected. 

Consider  the  case  where  Y  hes  h  equal  characteristic  roots 
which  are  much  larger  than  its  other  roots.  Then 


tr  ^  a  ^g1 


where  |_i  is  the  value  of  the  h  largest  characteristic  roots.  This 
results  in 


b  =  2  n 

Kg 

1 

°0  "  2  4g  rlh/2) 


and 


c.  =  0 
i 


i  >  0 
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Thus,  in  this  case  the  c_^  for  i  >  0  are  not  only  negligible,  but  ere 
identically  zero.  For  f(r  )  we  obtsin  then  the  chi-squared  distri¬ 
bution  with  h  degrees  of  freedom. 

If  two  unequal  characteristic  roots  dominate  the  others,  then  we 

have 


^3  =  -s(l  -  s)2 _ 

C0  (15/8)  +  2s  +  (17/4)  s2  +  2s3  +  (15/8)  s^ 


where  s  is  the  ratio  of  the  two  largest  roots.  This  function  has  a 
maximum  of  about  0.05  near  g  =  O.27.  Therefore,  for  our  purposes  we 
can  neglect  the  c_^  for  i  >0.  This  results  in 


f(r2)  =  1 


0  r(v2/02) 


V 

L  2 
0 


2  (v,/ 0)  -1 


exp 


vir 


{-7-} 


r2  >  0 


2 

which  is  of  the  form  of  the  chi-squared  distribution  with  2{v~Jo) 
degrees  of  freedom.  It  is  easy  to  show  that 

1  <  2(vl/0 )2  - 

tr  Y 


<  q. 
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V.  THE  CUMULATIVE  DISTRIBUTION  FUNCTION  OF  r 


The  cumulative  distribution  function  of  r  is  defined  by 


F  2(k2)  =  Prob.  (r  <  k)  =  Prob.  (r2  <  k2)  =(*  f(r2)  d(r2) 

r 


Using  the  above  result  for  f(r  )  this  becomes 


„  „  vl  *  r  V  l(Vo)  -1 
F  2<k  >  =  2  .  '  g,  g."  J0  —  exp 

r  a  r(v1  /0  )  u  L  o  J 


This  is  of  the  form  of  the  incomplete  gamma-function  which  has  been 

(4)  (5) 

tabulated  by  Pearson  and  Kendall  and  is  defined  as 


-  7Tp  ru  4 


-v  p  , 
e  v  dv 


P  :>  “  1 

u  >>  0 


Therefore,  using  this  notation  we  can  write 


2  Tk2  V1  2 

p  „(k2)  - 1 L#- ,  (— )  - 1 J 


where 


=  tr  Y  =  mean  squared  miss  distance 

2  2 

a  =  2  tr  Y 


F  2(k  )  is  plotted  in  Figs.  1  thru  4.  For  the  convenience  of  the 
r 

reader  the  same  function  is  shown  on  several  different  scales. 


If  tne  c. ;  i  >  0,  are  not  neglected,  then  clearly  F  2  can  be 

r 

represented  exactly  as  a  sum  of  incomplete  gamma-functions. 


Probability  (r<  k  )=  Fr2  (k2) 


1.0 


Probability  (r  <  k)  =  F r 2  ( K2 ) 
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k//o~ 

Fig.2  —  Cumulative  probability  distribution  of  r; 
Rayleigh  versus  linear  scale 


Probability  (r<k)  =  Fr2  (k 


k  X/cr 

Fig.3 — Cumulative  probability  distribution  of  r;  log-Rayleigh  versus  log  scale 


Probability  (r  <  Wd)  -  Fr2(k  ) 


0.9999 


0.999 


kZ/<r 

Fig.  4  —  Cumulative  probability  distribution  of  r2*  log-Rayleigh  versus  log  scale 
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VI.  APPLICATION  TO  TEE  PROBLEM  AT  HAND 

2 

In  order  to  determine  end  a  we  must  know  tr  Y  and  tr  Y  where 

Y  =  DMDT  +  Z 

Since  Y  is  at  most  a  3  x  3  matrix,  this  is  a  simple  and  straight¬ 
forward  calculation  even  in  the  most  general  case.  However,  many  times 
two  simplifications  may  be  made,  and  we  will  pursue  these. 

If  the  measurements  made  on  the  ballistic  missile  are  uncorrelated, 
then  K  is  diagonal  and  can  be  written 


Here  M  is  the  mean  squared  value  of  the  error  in  the  i-th  measurement. 
There  are  n  measurements  in  all.  If  further,  the  noise  in  the  AMM 
guidance  is  such  that  the  uncertainty  in  the  position  of  the  AMM 
at  time  t  has  spherical  symmetry  (e.g.,  when  only  the  rms  value  of 
the  AMM  miss  distance  is  known),  then  we  can  write 

z  =  (C/a)  I 

where  £  is  the  mean  squared  miss  distance  due  to  AMM  guidance  noise. 

* 

Using  these  substitutions  and  the  properties  of  the  trace  operator 
we  obta in 

tr  Y  =  Q  +  tr  DTDM 

and 

tr  Y2  =  C2/l  +  (2£/q)  tr  DTDM  +  tr  j^(DTDM)2] 
tr  ABC  =  tr  CAB  end  tr  (A  +  B)  =  tr  A  +  tr  B. 


lb 


These  can  "be  further  simplified  to 

n 

tr  Y  =  C  +  Y  \  \ 


and 


where 


i=l 


n 


n  n 


tr  Y2  =  QSh  *  (2 C/9)  l  A.  M.  +  l  l  B  M±: 


M. 

0  1  j 


“(  I  Dki 


i=l 


n,  .  =  b.. 


i=l  0-1 


hi  ko 


Ji 


,k=l 


Ai»I  Vs--/ 
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0=1 


and  D . .  is  the  j,i-th  element  of  the  sensitivity  coefficient  matrix, 
0  ^ 


given  by 


bw . 

D . .  =  -s-p- 

0i  am 


Here  w.  is  the  j-th  position-estimate  rectangular  coordinate  and  mu 

J 

is  the  i-th  measurement.  If  we  now  arrange  the  into  a  column  matrix 
/>rt=  the  into  a  column  matrix  A  =  jA^j-,  and  the  into  an 

=  Kj}» 


n  x  n  matrix  B 


then  we  can  write 


tr  Y  =  q£  +  A  /m 


and 

2  2  T  T 

tr  Y  =  q C  +  2 Q  A  M  +  /»TB/yn 

Since  A  and  B  are  easily  deteimined  from  the  sensitivity  coefficients, 
once  these  coefficients  are  determined  it  is  easy  to  calculate  =  tr  Y 
and  a  =  (2  tr  y2)1/2.  Then  any  of  Figs.  1  thru  4  can  he  used  to  obtain 

the  probability  that  r  <  k  as  a  function  of  k. 


VII.  SYSTEM  OPTIMI ZATI ON 


Often  the  geometry  and  the  kind  of  measurements  to  he  made  are 
difficult,  if  not  impossible,  for  a  system  designer  to  change.  When 
this  is  the  case  the  sensitivity  coefficients  and  the  resulting  A  and 
B  matrices  are  fixed.  The  designer  must  then  decide  how  much  money, 
weight,  and  equipment  to  allot  to  making  the  various  measurements  ac¬ 
curately. 

From  the  figures  it  is  evident  that  the  most  important  factor  in 

increasing  the  probability  that  the  AMM  conies  within  a  given  distance 

of  the  BM  is  the  reduction  of  o.  Accordingly,  an  important  quantity 

is  the  sensitivity  of  a  to  changes  in  the  M..  Since  it  is  usually  the 

J 

relative  size  of  fa  and  relative  sizes  of  the  RMS  measurement  errors 
which  are  of  interest,  we  will  let 


s  =  a(log  fc)  =  x  _  Jj: _ 

1  d(log/M±)  a  ^i  tr  Y2 


d(tr  Y2) 


dM. 


Now 


,  if  the  Si  are  arranged  in  a  column  matrix  S  =  -^S^j,  we  have 


S  = 


tr  Y 


MA  +  2  MBrtTl 


'] 


Denoting  by  A  /a  the  change  in  fa  caused  by  a  change  of  in  , 

we  can  write^  using  a  linear  extrapolation, 


hfc 

fo  : 


A/M, 


S. 


1  M 


Thus,  once  the  S  matrix  is  determined,  the  relative  importance  of  change 
in  the  accuracies  of  the  various  measurements  will  be  evident.  The  de¬ 
signer  can  then  proceed  in  a  systematic  way  to  make  changes  which  will 
optimize  the  system. 
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VIII.  EXAMPLE 


For  8n  example  we  consider  the  set  of  measurements  given  in  Table  1. 
For  a  particular  geometry  the  sensitivity  coefficient  matrix  for  these 


Table  1 


MEASUREMENTS  AM)  THEIR  ACCURACIES 


Measurement 

RMS  Error 

M. 

1 

Range,  R 

100  ft 

10,000  ft2 

Azimuth  angle ,  A 

50  4  rad 

2,500  (4  rad)2 

Elevation  angle,  E 

50  4  rad 

2,500  (4  rad)2 

Range  rate,  R 

10  ft/ sec 

100  (ft/ sec) 2 

Azimuth  angle  rate,  A 

5  4  red/sec 

25  (4  rad/sec)2 

Elevation  angle  rate,  E 

5  fj.  rad/ sec 

25  (4  rad/sec)2 

measurements  has  been  calculated  in  Ref.  1  (Case  IB,  p.  82)  to  be 


ft/ft 

ft/4  rad 

ft/4  rad 

ft/ ft/ flee 

ft/4  rad/aec 

ft/4  rad/ sec 

1.37 

-2.91 

-3.06 

2.78  x  102 

6.87  x  102 

-1.26  x  102 

0  » 

I.63  x  10"1 

1.07 

2.8+ 

-2.88  x  102 

6.72  x  102 

-8.29  x  10 

2.21 

-3.36  x  10"7 

-1.91  x  lCf1 

2.01+  x  10 

1.1+9  x  102 

1.07  x  103 

B 

A 

E 

R 

A 

E 

Using  a  simple  digital- computer  program  we  obtain  from  this 

— 

8.79 

1.20  x  10 
1.75  x  10 
l.bl  x  105 
9-^b  x  105 

1.17  X  iob 


A 
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end 


4. 6l  z  10 

I.36  x  10 

1.72  x  10 

1.44  x 

105 

1.90  X  106 

4.73  X  106 

1.36  x  10 

1.43  x  102 

2.02  x  102 

1.82  X 

106 

5.51  X  105 

4.73  x  10k 

B  • 

1.72  x  10 

2.02  x  102 

3.05  x  ID2 

2.80  X 

106 

4.94  x  10^ 

2.32  x  103 

1.44  x  105 

1.82  x  106 

2.80  x  lO6* 

2.56  x 

1010 

2.40  x  105 

1.02  x  108 

1.90  x  106 

5.51  x  105 

4.94  x  104 

2.40  x 

10 5 

8.94  x  1011 

2.49  x  10® 

4.73  x  106 

L 

4.73  X  10 

2.32  x  103 

1.02  x 

108 

2.49  x  108 

1.36  x  1012 

"We  will 

assume  the 

miss  distance  due 

to 

AMI''!  guidance 

noise  to 

have  e.  spherical  distribution  end  will  consider  three  different  values 

for  its  rms  value:  0,  1000  ft,  and  10,000  ft.  Thus,  we  have  £  =  0, 

6  8 

10  ,  and  10  .  Tnis  results  in 


6.91 

X 

107 

ft2 

c  = 

0 

tr  Y  =  < 

7-01 

X 

107 

ft2 

c  = 

106 

ft 

r 

H 

ON 

VO 

X 

108 

ft2 

c  = 

108 

ft 

r 

1.68 

X 

1015 

ft7 

c  = 

0 

2 

tr  Y  =  < 

1.72 

X 

1015 

h 

ft 

c  - 

106 

ft 

9.61 

X 

1015 

ft7 

c  = 

108 

ft1 

2(v1/a)2 


"2.84 
<  2.85 
^2-97 


/ff  = 


7610  ft 

<  7662  ft 

11,780  ft 


C  =  0 

6  2 

C  =  10  ft 

8  2 

£  =  10  ft 


C  =  0 

C  =  106  ft 
C  =  108  ft 


2 

2 
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C  =  0 

C  =  106 

1.08  x  10'3 

1.08  x  lo'3 

-4 

2.95  ^  10 

_k 

2.98  X  10 

_k 

4.21  x  10 

> 

-4 

4.27  x  10 

1.55  x  10_1 

1.57  x  10"1 

3.34  x  10_1 

3.34  x  10_1 

5.10  x  10_1 

5.07  x  10_1 

7  _k 
6.59  x  10 

R 

-4 

2.59  x  10 

A 

.  -4 

3.76  x  10 

TP 

U 

1.38  x  lo"1 

R 

-1 

•  . 

2.22  x  10 

A 

-1 

♦ 

2.91  x  10 

E 

From  the  curves  we  can  see  now  that  for  these  values  of ^/a  the  prob¬ 
ability  is  0.95  that  the  miss  distance  will  be  less  than 


13,500  ft  for  ■Jl 

13,500  ft  for  JZ  =  1,000  ft, 

end  20,600  ft  for  =  10,000  ft. 

From  the  S  matrices  we  note  that  is  about  100  times  as  sensi¬ 
tive  to  small  changes  in  the  accuracies  of  the  rate  measurements  as  it 
is  to  changes  in  the  accuracies  of  the  other  measurements.  According¬ 
ly,  we  consider  making  the  rate  measurements  twice  as  accurately  and 
making  the  others  only  half  as  accurately  (see  Table  2).  This  might 
result  in  a  system  of  about  the  same'cost."  This  change  should  reduce 
*/a  by  about 


.5  (1.55  +  3.34  +  5.10)  x  10"1  =  50$  for  £  =  0 

.5  (1.38  +  2.22  +  2.91)  x  10_1  =  32$  for  £  =  108 


end 
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Table  2 

MEASUREMENTS  AND  THEIR  REVISED  ACCURACIES 
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and  S  = 


3857  ft 
3959  ft 
9810  ft 


C  =  0 

C  =  106  ft2 
C  =  108  ft2 


C  =  0 


C  =  10 


C  =  10" 


1.71  x  lo-2 

1.69  x  10'2 

4.32  X  10"3 

4.80  x  lo" ^ 

4.98  x  10_3 

1.84  x  10" 3 

6.87  x  10“3 

y 

7.14  x  10~3 

2.68  x  lo“3 

I.58  x  10"1 

1.64  x  10_1 

6.16  x  10’ 2 

3.21  x  10'1 

3.21  x  lo"1 

9-28  x  10‘2 

4.93  x  10"1 

4.84  x  lo”1 

1.17  x  Kf1 

R 

A 


E 


R 

A 

E 


The  new  values  of  ^/o  check  fairly  well  with  the  estimates  obtained 
above.  Tne  correspondence  is  not  perfect  because  the  S  matrix  is  not 
strictly  independent  of  the  M  (i.e.,  the  linear  extrapolation  is  not 
perfect).  Now  the  probability  is  0.95  that  the  miss  distance  will  be 
less  then 

6840  ft  for  ^  -  0, 

7010  ft  for  Jq  =  1000  ft, 
and  17,400  ft  for  =  10,000  ft. 

This  procedure  of  observing  S  and  then  making  changes  which  will 
improve  the  performance  while  maintaining  or  reducing  the  cost  of  the 
system  can  be  continued  until  a  satisfactory  design  is  obtained. 
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